The pulling force of a single DNA molecule condensed by spermidine 
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In a recent experiment, a single DNA double helix is stretched and relaxed in the presence of 
spermidine, a short positive polyelectrolyte, and the pulling force is measured as a function of DNA 
extension. In a certain range of spermidine concentration, a force plateau appears whose value 
shows maximum as a function of spermidine concentration. We present a quantitative theory of this 
plateau force based on the theory of reentrant condensation and derive almost parabolic behavior 
of the plateau force as a function of the logarithm of the spermidine concentration in the range of 
condensation. Our result is in good agreement with experimental data. 

PACS numbers: 87.15.La, 61.41.+e, 64.70.-p, 87.15.He 
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DNA condensation by strongly positively charged pro- 
teins, histones, is used by nature for compaction of DNA 
in cell nucleus (DNA is strongly negatively charged). 
Positive proteins, protamines, are used for additional 
compaction of DNA in the sperm . Gene therapy uses 
complexes of DNA with long positive polyelectrolytes or 
other macrocations. The net charge of these complexes 
can be positive and therefore they are not repelled by 
negative cell membrane in the course of gene delivery. 
Thus, condensation of DNA by positive macrocations is 
an extremely important physical phenomenon. It is in- 
tensively studied in the model system of double helix 
DNA with spermidine, a flexible polymer with length 15 
A and charge +3. It is known that in a dilute DNA 
solution at some concentration of spermidine, s = s^, 
each long DNA molecule self- condenses into a toroid [2j • 
When s grows farther, at a much larger concentration, 
s = Sd, DNA dissolves back This phenomenon is 
called reentrant condensation and has got theoretical ex- 
planations 




FIG. 1: Schematic illustration of the experimental setup. 
Partly condensed long DNA is pulled by two beads. 

Recently, a new single molecule technique has been 
used to measure the force necessary to pull a single 
DNA double helix from a toroidal condensate of DNA 6] 
(Fig. fT)l . In the experiment, extremely small concentra- 
tion of DNA (2 nM of nucleotides) is dissolved in wa- 
ter with certain concentrations of spermidine (varying 



from 200 fiM to 200 mM). Dual-trap tweezers are used 
to stretch a single DNA molecule tethered between two 
protein-coated polystyrene beads. DNA is first stretched 
and then relaxed. During these processes, the pulling 
force / is measured. Experiments are repeated at vari- 
ous spermidine concentrations. It is observed that in the 
interval s c < s < s<j, there is a fairly large range of DNA 
extension, x, where force / is significantly larger than the 
force described in the wormlike chain model Q. More- 
over, in this range of x, f is almost a constant (a force 
plateau, see Fig. |2J). The plateau value of the force, / p , 
is plotted as a function of s in Fig. [3] 







FIG. 2: The pulling force / measured as a function of exten- 
sion x is schematically shown by the solid line for s c < s < Sd, 
when DNA is condensed by spermidine. At x < xo, the force 
f p is a constant whose value depends on s. The dashed line 
shows the wormlike chain behavior of / in the absence of DNA 
condensation (s < s c or s > s^). At x > xo, the condensate is 
completely eliminated by the pulling force and the wormlikc 
chain behavior is recovered. 

In this paper, we suggest a quantitative theory of the 
/ p (s) curve based on the theory of reentrant condensa- 
tion [j| . We assume that the whole process is slow enough 
so that the system is always in equilibrium. Then in zero 
order approximation, the force needed to detach DNA 
from the condensate is just the free energy difference per 
unit length between the free coil state and the condensed 
state of the DNA-spermidine complex. We call it /q. 
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The pulling force f p measured in the experiment is some- 
what larger than /o since the detached part of the DNA 
molecule is stretched, and therefore loses entropy with 
respect of its free coil state. In the end of the paper, 
we show that this effect adds a small correction to the 
pulling force, 

when ksT/lfo <C 1. Here I is the persistence length of 
DNA. In this paper, we first calculate fo and then f p us- 
ing Eq. JJJ . Our result is shown in FigEJ It demonstrates 
reasonably good agreement with experimental data. 
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FIG. 3: Comparison between the theoretical result (the solid 
line) and experimental data (points with error bars taken from 
Fig. 2 of Ref. 0) for the pulling force f p as a function of 
spermidine concentration s. fo is shown by the dashed line. 

Let us start from the theory of reentrant condensa- 
tion of DNA with spermidine without any external forces. 
In the solution, positively charged spermidine molecules 
are absorbed on the negatively charged surface of dou- 
ble helix DNA. Since spermidine molecules are strongly 
charged (+3), they form a two dimensional strongly cor- 
related liquid on the surface of the DNA molecule. When 
a new spermidine molecule approaches the DNA helix, it 
repels already absorbed spermidine molecules and creates 
an electrostatic image of itself, similar to the image on 
a conventional metallic surface. Attraction to the image 
leads to an additional negative chemical potential fj, c for 
spermidine molecules on the surface of DNA. As a result, 
if s exceeds some concentration, sq, charge of the DNA- 
spermidine complex changes sign and becomes positive, 
i.e., spermidine molecules overcharge DNA (sec review of 
the theory of charge inversion in Ref. || ) . 

Adsorbed layers of the correlated spermidine liquid 
also lead to self-attraction of DNA. In the spot where two 
turns of DNA touch each other, the surface density of the 
correlated liquid is doubled and the correlation energy 



per spermidine molecule is reduced [9|- This short range 
force of self-attraction leads to condensation of DNA if 
the macroscopic Coulomb repulsion is not very strong. 
As a result, in the vicinity of the neutralizing concen- 
tration so ( s c < so < Sd), DNA-spermidine complex is 
condensed into a toroid. 

The condition of equilibrium between spermidine 
molecules in the bulk and on the surface of DNA has 
a form 

(j, c + Ze<j> = k B Tln(svo), (2) 

where Z = +3 is the valence of spermidine, <f> is the 
electric potential on the surface of the DNA-spermidine 
complex, and vo is the normalizing volume of spermidine. 
The left hand side is the chemical potential of a spermi- 
dine molecule in the complex, the right hand side is its 
chemical potential in the bulk of the solution. Notice that 
the entropy of spermidine can be expressed through the 
total concentration s because absorbtion of spermidine 
by DNA practically does not change s (the concentration 
of DNA nucleotides is at least 10 5 times smaller than s). 
In the vicinity of sq (s c < s < Sd), the net charge of 
the complex is small and the spermidine concentration 
on the DNA surface does not change much. Since \i c is 
determined by this concentration, it is approximately a 
constant. According to Eq. J5J, when <j> — 0, i.e., the 
complex is neutral, there is a simple relation between /j, c 
and so, 

S0 = fexp(-JMY (3 ) 
v \ k B T J 

Now we calculate the free energy of the charged com- 
plex in its coil state. For this purpose, we treat the DNA- 
spermidine complex as a capacitor with the capacitance 
per unit length C . Suppose s is such that the complex is 
overcharged (so < s < s<j). Overcharging enhances the 
Coulomb self-energy and reduces both the correlation en- 
ergy and the entropy of spermidine molecules. Taking the 
free energy of a neutral coil as zero, and using Eq. J2J, 
we get the free energy of the complex per unit length 

h = + §7 Mc ~ ^ B THm) = -\c<P, (4) 

where the three terms allow for the three parts of the 
free energy mentioned above, and C<p/Ze is the number 
of spermidine molecules overcharging the complex. The 
final expression is the same as the free energy of a capaci- 
tor kept under a constant voltage <j>. It is easy to see that 
this expression is also true for an undercharged complex. 
Since the DNA-spermidine complex can be considered as 
a long cylinder, C is given by expression 

c D ( 51 

21n(l + r s /i?)' 1 ' 
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where D = 80 is the dielectric constant of water, r s is 
the Debye-Hiickel screening radius, and R = 10 A is the 
radius of the double helix cross section. 

We first assume that r s R- In this case, the complex 
in the condensed state is practically neutral. Indeed, if 
the condensed complex were charged, the Coulomb-self 
energy of the macroscopic condensate would be too large 
to hold it. We define phenomenological parameter e < 
as the free energy per unit length of the condensed com- 
plex calculated from the free energy of a neutral coil. It 
includes the gain of correlation energy in the spots where 
two turns of DNA touch each other, and also the loss of 
certain entropic elasticity of the coil in the condensate. 

When the complex goes from the neutral condensed 
state to the charged coil state, the free energy increment 
per unit length is 



fo goes to zero, 



h = -\c4> 2 -e. 



(6) 
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FIG. 4: The zero order approximation to the pulling force, 
fo , as a function of In s given by an inverted parabola. The 
maximum value of fo, \e\, is achieved at s = so where the free 
DNA-spermidine complex is neutral. s c> d are the concentra- 
tions of spermidine at which the condensate dissolves. 

If r s is comparable with R (as in experiment 6J), the 
approximation used above that the condensate is neutral 
has to be revised and the Coulomb energy of the conden- 
sate should be taken into account. If we introduce C as 
the effective capacitance per unit length of the complex 
in the condensate, similarly to Eq. the free energy 
per unit length is just — C'(j) 2 /2. Accordingly, in Eq. 10, 
C should be replaced by C - C . Using Eqs. @ and 0, 
we rewrite Eq. as 



fo = \e\ - 



(C - C')k 2 B T 2 
2Z 2 e 2 



In*-. 

so 



(7) 



The function /o(lns) is shown in Fig. 0] by an inverted 
parabola. According to Eq. JJJ, the maximum value fo = 
\e\ is achieved at s — sq where the free complex is neutral. 
Also, we get the two concentrations s c and Sd at which 



s exp =F- 



2\e\Z 2 e 2 



kaT V C-C 



(8) 



where the upper (lower) sign corresponds to the first (sec- 
ond) subscript of s c .d- Exactly at these two concentra- 
tions, one can see transition from the coil state to the 
condensed state or vice versa in light scattering experi- 
ments with solutions of DNA and spermidine |5j • 

In order to calculate C", we assume that the condensate 
is macroscopic and densely packed. Then 



pe r/rs A 2 , 4:Trr 2 p 
-A-Kr 2 dr = 



Dr 



D 



(9) 



where p is the charge density of the condensate. The 
charge of the complex per unit length is irR 2 p/a where 
a = 0.91 is the filling factor for the hexagonal dense 
packing of cylinders. This gives 



c _ p *R 2 

4> a 
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We see that when r s 3> R, we can drop C' in C — C 
because C <^L C, i.e., the condensate is almost neutral. 

As we mentioned before, fo is only the zero order ap- 
proximation to the pulling force. The pulling force f p is 
given by Eq. Q • We need four experimental parameters 
to calculate it. The spermidine concentrations s c .d at 
which the condensate dissolves, the average DNA persis- 
tence length I, and the average Debye-Hiickel screening 
radius r s . Following Ref. 0], we take s c = 0.35 mM, 
s d = 150 mM. Therefore s = 7.2 mM (Using Eq. ©). 
We use I = 500 A corresponding to the persistence length 
of a neutral DNA coil. At so = 7.2 mM, spermidine and 
its counter ion contribute to r s even more than monova- 
lent salt (10 mM). Treating every spermidine molecule as 
a point-like trivalent ion, we g et r, s = 13 A from the stan- 
dard Debye-Hiickel expression [lfj . Calculating C and C 
according to Eqs. JSJ and (|10|l . we finally get |e| = 0.11 
k B T/bp (lbp = 3.4 1) from Eq. ©. 

For these parameters, fo and f p are shown together 
with the experimental data in Fig. [31 We see that our 
result for f p agrees pretty well with the experimental 
data. Notice that correction f p — fo is much smaller than 
fo in the range of experimental data. This justifies the 
use of the perturbative result Eq. 

Finally, let us derive Eq. For this purpose, we 

consider the detached part of DNA as a wormlike chain 
with the contour length L and end-to-end distance x 
(see Fig. nj, while the contour length of the whole DNA 
molecule is L - We also assume that the system is always 
in equilibrium during the whole process. The free energy 
of the system is 

F(x,L) = -{L -L)fo + F e 
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-(L - L)f 



L 



2L 4(1 -x/L) 



Solving Eq. (|13|l in the same limit, we get 



(11) 



where — (Lq — L)fo is the free energy of the condensed 
part of DNA, and F e is related to the entropic elas- 
ticity of the detached part of DNA. The expression 
for F e is obtained for a free wormlike chain with fixed 
contour length L and end-to-end distance x: F e (x) = 
Jq f e {x')dx' + const.. Here f e is the force needed to keep 
the end-to-end distance x for a free wormlike chain with 
length L and is given by an interpolation formula 



fe 



k B T 
I 



1 



x 1 
I ~ 4 + 4(1 -x/L) 2 



(12) 



The energy zero point has been chosen at x = 0, L = L , 
i.e., at the free energy of a DNA-spermidine complex in 
its coil state. 

For a given x, the contour length L of detached DNA 
can be found from the condition of minimal free energy, 
namely, dF/dL\ x = 0. If we define a = x/L < 1, this 
minimum condition can be written as 



2a 2 (l -a) 



a 2 - 



4Z/o 



(1-a) 2 = 0. 



(13) 



Finding a from this equation and putting L = x/a back 
to Eq. 1|11|) . we calculate the magnitude of the pulling 

force 



dF 

dx 



fo 
a 



k B T 
I 



1 



4 4(1 -a) 



(14) 



Notice that a does not depend on x and L (see 
Eq. §l'S\i ). Therefore, according to the definition x = aL, 
the two lengths x and L are always proportional to each 
other. This result was also obtained in Ref. [llj using 
a different approach. From Eq. (|14|l it is easy to see 
that f p > fo since a < 1. There are two contributions 
to f p — fo. The first is a geometrical effect (the term 
fo/a in Eq. I|14|l % ). Namely, to change the end-to-end dis- 
tance by Ax, a larger length AL — Ax /a must be pulled 
out of the condensate. The second contribution is the 
free energy increment of the detached part of DNA, F e . 
These two contributions are equal in the limiting case of 
k B T < Ifo (see Eqs. (15}, (JTHJ) and (J)). Since f p does 
not depend on x, we get a force plateau with increasing x, 
as observed in the experiment (see Fig.|2J. This plateau 
ends when x reaches xo = aLo and all DNA is pulled out 
of the condensate. At x > xq, the pulling force starts to 
increase with x according to Eq. I|12l) 12j . 

When k B T/lf < 1, we have 1 - a < 1, and Eq. Q 
becomes 



, _ fo 

h ~ a '• 41(1 



k B T 



■a) 



(15) 




(16) 



Substituting Eq. (|16|) in Eq. I|15l) . we arrive at our fi- 
nal result Eq. QJ. It gives the pulling force f p in the 
first order perturbation theory in the small parameter 
\/k B T /Ifo. In experiment |6j, we have k B T/lfo = 0.06 
at s — so. Therefore, Eq. can be used in almost all 
the range s c < s < s^. At s close to s c .d, where fo — > 0, 
the perturbation theory fails. In principle, we can calcu- 
late f p numerically using Eqs. (|13|l and (|14fl in the whole 
interval s c < s < s<j. 
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Strictly speaking, at x —> xo, our theory is not valid 
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can not be ignored in this limit). Therefore at x close 
to xo, the pulling force / first decreases a little bit with 
increasing x, and then goes back to the wormlike chain 
behaivor. 



